Non-equilibrium steady state in a periodically driven Kondo model 
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We investigate the Kondo model with time-dependent couplings that are periodically switched on 
and off. On the Toulouse line we derive exact analytical results for the spin dynamics in the steady 
state that builds up after an infinite number of switching periods. Remarkably, the universal long- 
time behavior of the spin-spin correlation function remains completely unaffected by the driving. 
In the limit of slow driving the dynamics becomes equivalent to that of a single interaction quench. 
In the limit of fast driving it is shown that the steady state cannot be described by some effective 
equilibrium Hamiltonian due to the observation that an incautious implementation of the Trotter 
formula is not correct. As a consequence, the steady state in the limit of fast switching serves as an 
example for the emergence of new quantum states not accessible in equilibrium. 
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I. INTRODUCTION 

Recent progress in experiments stimulated the inter- 
est in non-equilibrium phenomena of interacting many- 
particle systems. Cold atoms trapped in optical lattices 
offer the possibility of studying the time evolution of 
quantum many-body systems with time-dependent sys- 
tem parameters.^ Due to the excellent isolation from the 
environment the non-equilibrium dynamics of these sys- 
tems are accessible with negligible decoherence over long 
times. 

Nanostructures such as quantum dots provide 
the framework to examine experimentally the non- 
equilibrium dynamics in quantum impurity models. Most 
importantly for the present work, quantum dots can act 
as magnetic impurities displaying Kondo physicsP In 
contrast to impurities in a bulk sample, unscreened elec- 
trical or magnetic fields can be applied directly such that 
the local system parameters can be varied in time by 
choosing appropriate time-dependent fieldsP 

The possibility to experimentally study the proper- 
ties of interacting many-body systems out of equilibrium 
motivated numerous analytical and numerical theoretical 
treatments. Most of the activities have been concentrat- 
ing on interaction quenches in various model systemsP^ 
For periodically driven interacting many-body systems, 
however, less results are known. Recently, periodic time- 
dependent Falicov-Kimball models have been investi- 
gated in the limit of infinite dimensions by using dy- 
namical mean-field theoryP Considerable activity in the 
field of time-dependent quantum impurity models led 
to a nu mber of works on periodically driven Anderson 
impuritjEEHl an d Kondo modelg^Sl. 

As the Kondo model is the paradigm model for 
strongly correlated systems, it is of particular interest 
in the field of non-equilibrium phenomena. In equilib- 
rium, the Kondo effect emerges from the interaction of a 
localized spin degree of freedom with a bath of surround- 
ing electrons. At sufficiently low temperatures, this bath 
of itinerant electrons develops a localized spin polariza- 



tion cloud in the vicinity of the local spin, the so-called 
Kondo cloud, providing a mechanism to screen the local 
magnetic moment. In the zero temperature limit, the 
screening becomes dominant leading to the emergence 
of a bound state. The surrounding spin polarization 
cloud is tied to the local spin establishing the so-called 
Kondo singlet with an associated binding energy Tk, the 
Kondo temperature. The Kondo effect manifests itself 
most prominently in the Kondo resonance, a sharp peak 
in the local density of states that is pinned exactly at the 
Fermi energy. As the Kondo effect is a coherent many- 
body phenomenon, the question arises how it is affected 
in a non-equilibrium setting. 

Due to the complexity of many-body systems out of 
equilibrium, it is instructive to investigate those cases 
where exact nonperturbative solutions are accessible. In 
this context, it is of particular interest that the Kondo 
model as a paradigm model for strongly correlated sys- 
tems exhibits a special line in parameter space, the 
Toulouse limit, where it becomes exactly solvable.^ The 
Toulouse limit displays many generic and universal prop- 
erties of the strong coupling limit of the Kondo model in 
equilibrium as well as for interaction quenches^. The lo- 
cal spin dynamics, for example, that is also investigated 
in this work, is well described whereas other universal 
quantities such as the Wilson ratio explicitly depend on 
the anisotropy. 

The exact solvability of the Kondo model in the 
Toulouse limit is used in this work to investigate non- 
perturbatively a steady state that is generated by pe- 
riodically switching on and off the interaction at zero 
temperature. This steady state is characterized by an- 
alyzing exact analytical results for the local dynamical 
quantities, that is the magnetization of the impurity spin, 
the spin-spin correlation function and the dynamical spin 
susceptibility. 

A system that is driven by an external force approaches 
a steady state if the amount of energy that is provided 
to the system does not lead to an overheating, as it may 
happen for systems with an unbounded spectrum. Since 
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FIG. 1: Time dependence of the perpendicular J±(t) and 
parallel coupling Jy (t) in the anisotropic Kondo Hamiltonian. 
For times t < all couplings are zero whereas for t > they 
are switched on and off periodically with period r. 



these steady states emerge from a non-equilibrium set- 
ting, equilibrium thermodynamics is not applicable for 
their description. As a consequence, these states display 
new properties that are not accessible by exciting the sys- 
tem thermally. Their characterization, however, poses a 
new challenge. Recently, there have been attempts to 
assign effective thermodynamic quantities such as effec- 
tive temperatures to desc ribe th e properties of systems 
in a nonequilibrium setup EElZISI As will be shown in this 
paper, a characterization of the present steady state in 
terms of an effective temperature is not possible. The 
excitations that are created by the periodic driving are 
fundamentally different from those induced by temper- 
ature. A finite temperature smears the Fermi surface 
whereas the periodic driving leads to an excitation spec- 
trum of discrete character with excitations of multiples 
of the driving frequency. 

Another question arising in the context of driven sys- 
tems is whether the universality in equilibrium sys- 
tems such as the Kondo model extends to the non- 
equilibrium case, whether new universal properties 
emerge and whether equilibrium qua ntitie s remain mean- 
ingful. Kaminski and coworker ;] 11 ! 12 !, for example, 
proposed a universal description for the conductance 
through a Kondo impurity for a periodically driven 
Kondo model. In the present setting, the spin-spin cor- 
relation function displays a universal long-time behavior 
that is completely independent of the driving. This uni- 
versality originates from the fact that the low-energy ex- 
citations in the immediate vicinity of the Fermi level that 
are relevant for the long-time behavior are unaffected by 
the periodic driving due to the discrete character of the 
excitation spectrum. 

This paper is organized as follows. In Sec. II, the 
model Hamiltonian, a time-dependent Kondo model in 
the Toulouse limit, is introduced and mapped onto a 
quadratic effective Hamiltonian. The method used to de- 
termine the time evolution in the periodic driving setup 
is presented in Sec. III. The result for the magnetization 
of the impurity spin is shown in Sec. IV. Sec. V is devoted 
to a detailed analysis of the spin-spin correlation function 
and the results for the dynamical spin susceptibility are 
presented in Sec. VI. 



II. PERIODIC TIME-DEPENDENT KONDO 
MODEL 

Consider a local spin whose exchange interaction with 
the surrounding electrons is switched on and off period- 
ically, which generates a Kondo Hamiltonian with time- 
dependent couplings. For convenience we allow for an 
anisotropy in the exchange interaction leading to differ- 
ent couplings in z direction, J z = Jii , and in the xy-plane, 

J X Jy J_\_- 

H = J2 k:c\ a c ka . + ]T^E :^(0K-'^ Q ,(0):. 

ka i a, a' 

(i) 

The operator c^, Q creates an electron with wave vector 
k and spin a =f? i i n t he reservoir. The colons : . . . : 
denote normal ordering relative to the Fermi sea. The 
local spin operator S with components Si,i — x,y,z, 
is coupled to the local spin density of the conduction- 
band electrons whose components are determined by the 
Pauli matrices er^. The electron's dispersion relation has 
been linearized around the Fermi level and energies are 
measured in units of vp relative to the Fermi energy, i.e., 
vp = 1 and Ep = 0. As the local scatterer is assumed to 
be pointlike, only s-wave scattering occurs rendering the 
problem to be effectively one-dimensionalP^ 

For negative times, the system is prepared in one of the 
ground states l^o) of the noninteracting problem that 
are product states of the Fermi sea for the conduction- 
band electrons and a wave function for the local spin. As 
depicted in Fig. [T] at time t — 0, the periodic driving 
process starts by switching on the interaction. After half 
of the period r, t — t/2, the interaction is switched off 
until t = t. Afterwards, this procedure is continued, until 
after an infinite number of periods a steady state builds 
up, a state in which all real-time correlation functions are 
invariant under a discrete time shift of one period in all 
their time arguments. The energy scale associated with 
the periodic switching is the driving frequency Q = 2it/t. 

Due to the time-dependence of the Hamiltonian, en- 
ergy is not a conserved quantity. Moreover, each quench 
that is performed excites the system such that one can 
expect that after an infinite number of switchings, in the 
steady state, an infinite amount of energy is pumped into 
the system. This is indeed the case in the present setup. 
Therefore, a dissipation mechanism is needed that pre- 
vents the system from overheating. As has been shown by 
Doyon and AndreP^, the conduction band in the Kondo 
model itself can serve as a bath if and only if it is taken 
as infinitely large. Therefore, a definite order of taking 
limits has to be prescribed, namely, the thermodynamic 
limit has to be taken before the limit of long times: 

lim lim . (2) 

t— >oo L— too 

In this way, as expected the conduction band is prevented 
from overheating induced by a single impurity. 
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One possible way of experimentally realizing the pe- 
riodic switch on and off of the Kondo interaction in 
a quantum dot is the following: consider a quantum 
dot with a local single-particle level at energy Ed and 
a large on-site interaction U in the Kondo regime where 
£d < EF,£d + U >£f and |e<j — £f|, \e<i + U — £f\ > T, T. 
Here, T denotes the broadening of the local level Ed in the 
quantum dot and T the temperature. Via the Schrieffer- 
Wolff transformation, the corresponding Anderson im- 
purity model can be mapped onto a Kondo model with 

" P Fol- 



identity reduces to: 



an exchange coupling J oc T 

lowing a suggestion by Nordlander and coworker^, con- 
sider the case where the local single-particle level Ed al- 
ternates between two different values Edi and Ed2 with 
|ed2 — £f\ ^ \sdi — sf\ where for each Ed\ and Ed2 the 
quantum dot is assumed to be in the Kondo regime. 
Then, the Kondo exchange coupling Ji corresponding to 
Ed2 is much smaller than J\ . The associated Kondo tem- 
perature fcsTfj-2 = Dy/ poJ2e~ 1 /( J ' 2Po \ po is the density of 
state at the Fermi level and D a high-energy cutoff, van- 
ishes exponentially such that J2 can be set equal to zero. 
As a result, the corresponding Kondo model becomes 
time-dependent with an exchange interaction J\ that is 
switched on and off periodically. This can be shown rig- 
orously by performing a time-dependent Schrieffer- Wolff 
transformation!^^ 

The periodic driving in a quantum dot sets an up- 
per bound on the driving frequency ft. As the Kondo 
Hamiltonian requires strict single occupancy, the driv- 
ing has to be small enough not to induce charge fluctu- 
ations on the dot caused by hopping processes between 
the central region and the conduction band by absorb- 
ing or emitting quanta of the driving frequency, i.e., 

n^\E d \,Ed + um^ 

At zero temperature in equilibrium, the screening of 
the local magnetic moment by the conduction-band elec- 
trons becomes dominant leading to the emergence of 
a bound state called the Kondo singlet. Many uni- 
versal features of the equilibrium Kondo model in this 
strong coupling limit arc well described by the anisotropic 
Kondo Hamiltonian in the Toulouse limit that corre- 
sponds to a special line in parameter space of the Hamil- 
tonian in Eq. |lj where Jy = 2 — y/2. For this value of the 
parallel coupling, the Hamiltonian can be mapped onto 
an exactly solvable quadratic noninteracting resonant- 
level model using bosonization and refermionizationP^ 
Recently, it was shown by Lobaskin and KehreirP that 
these methods can also be adopted to an exact solution 
of an interaction quench scenario. The only difference 
is an additional potential scattering term in the effective 
Hamiltonian. 



(3) 



where the bosonic field <fr a (x) — 
aq / 2 / ^/riq is related to 



e lqx b\ a ]e 



the fermionic densi- 
ties p a (x) = -.^(x^a^x): = d x (j) a (x) and a^ 1 is 
an ultraviolet cutoff. The bosonic operators b qa — 

^1 \[™q Sfe c \+ qa c ku create a superposition of particle- 
hole pairs with momentum transfer q = 2nn q /L > 0. 
Here, L denotes the system size. The Klein factor F a ac- 
counts for the annihilation of one electron as this cannot 
be achieved by the bosonic field <j> a (x). 

By performing a sequence of unitary transformations, 
the Kondo Hamiltonian in the Toulouse limit can be 
simplified tremendously. First, the spin and charge 
(s,c) degrees of freedom are separated by defining the 
bosonic fields (j) s {x) = [4>t( x ) — ^(aOJ/v^ and (f> c (x) — 
[<l>t(x) + (f)±(x)]/y2. The charge sector of the anisotropic 
Kondo Hamiltonian is decoupled from the local spin and 
reduces to a collection of uncoupled harmonic oscilla- 
tors. Therefore it will be omitted from now on. In 
the first half period, the interaction part in the spin sec- 
tor J {l /y2d x <j> s (x)S z + Jj{2a)[F\F i ^^S- + h.c] 
is modified by an Emery-Kivelson transformation U = 
e n<Ps(x)s^ 7 = 1 /2_ i ) to J ± /(2a)[F\F i e^^S- + h.c] 

in the Toulouse limit where Jy = 2 — v2- For the 
second half period, the Emery-Kivelson transformation 
generates a scattering term oc (S z (t))d x ^> s (0), whose 
strength depends on the instantaneous magnetization of 
the impurity spin (S z (t)}W The exponentials appearing 
in the transformed interaction part can be refermion- 
ized by introducing new spinless fermionic fields ty(x) = 
a ~i/2p^ e -up B (x) usm g the inverse of the bosonization 

identity where F s = F^F-fW^ One may think of the 

*5>'(x) fields as creating spin excitations at point x in 
the reservoir. Moreover, another unitary transformation 
U2 = e l7rA ^ Si has to be imposed in order to arrive at a 
completely refermionized Hamiltonian and to ensure cor- 
rect anticommutation relations for all operatorspSl Here, 
Af s = g[iVf — N±] measures the total spin polarization 
of the conduction-band electrons. By defining the oper- 
ator d = e~ m ^ a ~ s ^S- and its Hermitian conjugate d) 
as well as by performing a mode expansion for the new 
fermionic fields c k — (2-kL)^ 1 / 2 J dx^{x)e lkx , one arrives 
at the following Hamiltonian: 



H= Y,k:clc k :+git)(S z {t))^: 



c\c k > : 



(4) 



The bosonization technique establishes a bosonic rep- 
resentation of fermionic fields ^ a (x) in one dimension 
called the bosonization identity, see RefP^ for a recent 
review. In the thermodynamic limit, the bosonization 



where g(t) = gQ{- sin(fit)), V(t) = F0(sin(Qi)), V = 
J ± ^Tr/2aL, g = (1 - V2)ir/L for t > 0, and g(t) = g, 
V(t) = for t < 0. For times Nt < t < Nt + r/2 
the Hamiltonian is a resonant-level model. For times 
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Nt + t/2 < t < Nt + t the dynamics are governed 
by a potential scattering Hamiltonian and the local d 
operators do not evolve in time such that (S z (t)) = 
(S z (Nt + t/2)). As it will turn out in the following 
analysis, the intermediate time evolution with the po- 
tential scattering Hamiltonian has no influence on the 
local spin dynamics such as the magnetization (S z (t)) at 
all. Therefore, it is not necessary to solve the dynamics 
in a self-consistent way. 

The Kondo scale can be connected to the parameters of 
the resonant-level model via the impurity contribution to 
the Sommerfeld coefficient in the specific heatP C imp = 
7 imp T where 7, mp = wtt 2 /3T k and w = 0.4128 is the 
Wilson number. In this way, the Kondo temperature Tk 
is determined by Tk = ttwA where A = V 2 L/2 is the 
hybridization function. 

The functional dependence between the spinless 
fermions c k and the conduction-band electrons c ka is 
highly nonlinear and nontrivial. The local spin observ- 
able S z , however, commutes with all unitary transforma- 
tions and can be connected to operators of the effective 
Hamiltonian in a simple way: 



S z = dU - 



(5) 



This relation allows to analytically calculate correlation 
functions that involve the S z observable such as the mag- 
netization of the impurity spin P{t) = (S z (t)), the spin- 
spin correlation function (S z (t)S z (t')} and the dynamical 
spin susceptibility x"(A e ): as it will be done in this work 
for the periodic driving setup. 



with a unitary matrix G. Despite the complexity of 
time evolution for time-dependent Hamiltonians, the pe- 
riodic driving as it is considered in this work involves 
two time slices, during which the Hamiltonian is con- 
stant. For each half period, the dynamics is determined 
by the Green's functions for a potential scattering Q^> 
and a resonant-level model Q. The Green's function for 
the potential scattering Hamiltonian varies for different 
periods, as the strength of the scatterer depends on the 
instantaneous magnetization of the impurity spin. There- 
fore, the Green's functions are labeled by an additional 
superscript where n stands for the number of the period 
after starting the periodic driving. The time evolution 
over one period transforms the single-particle operators 
in the following way: 



P) 

defining the unitary matrices M. (n) that are obtained by 
matrix multiplication. Note that the order of time evo- 
lution in the Hcisenberg picture is opposite to the order 
of time evolution in the Schrodinger picture: the opera- 
tors are first evolved according to the potential scattering 
Hamiltonian although it acts in the second half period. In 
this formulation, the problem of long-time evolution re- 
duces to a matrix multiplication problem since the long- 
time evolution is completely determined by evolving the 
single-particle operators over multiple periods: 



ci(Nt) 



M 



(AT) 



[M(1)...M(N)] U > 



III. TIME EVOLUTION 

As the Hamiltonian in Eq. Q is quadratic, the time 
evolution of the single-particle operators c k and d is 
entirely determined by the Green's functions Gw (t) — 
8(t)({ci(t),cU): 



(-1 



{t)=Y,G w (t)c Vl l,l' = k,d 



(6) 



.(8) 

Finding the matrix elements of an arbitrary multiplica- 
tion of some matrices can be a difficult task. The matri- 
ces M(n), however, display nice mathematical properties 
such that an analytical calculation can be carried out as 
is shown in Appendix |Xj The goal is to derive solvable 
recursion formulas by partially performing summations 
over intermediate indices. This yields the following rela- 
tions: 



M%>=e 



a -iNkr c -NAr/2 „-iNkr „-NAr/2 

-NAr/2 U AO _ , a e ~ e ' kA^-U 6 ~ 6 

) M dk - M dk —7- awo » M kd ~ Mkd- 



M{ N J = 5 kk ,e 



! N k 7 



&kk' 



e -iNkr _ e -iNk 



e " w — e 

-NAt/2 



p-ikr _ p-A-r/2 ' 
MkdMdk 

r + 



e 

^—iNkr 



,-At/2 



?—ikr ik'r 1 \& — ik 



'Tll e -ikT _ e -Ar/2j 



-iNk t 



\^—ikr g— At/2J ik'r g— At/2 f - . .. : _ ( — > r. . ( — > , _ . - 



-ik'r p— ikr] \p—ik'r p — At/21 



(AT) 
kk' ■ 



(9) 



The precise definition of the functions appearing in these ments M.\y can be interpreted as the probability am- 
relations can be found in Appendix [A] The matrix ele- pfitudes for a fermion in a single-particle state \l) = c\\) 
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FIG. 2: Probability l-M^?" 1 ! 2 for a local d-fermion to decay 
into a bath state k after an infinite number of periods. As 
a reference, the equilibrium curve for a resonant-level model 
with the same Kondo temperature Tk is included. The pa- 
rameter 77 = tTk = tI%k compares the speed of switching r 
with the internal time scale tn = 1/Tk- 



to transform into after N periods r. Here, |) denotes 
the true vacuum without any fermion. In Fig. [2j plots 
for |A4^'| 2 are shown, that is the probability for a d 
fermion to decay into a k fermion after an infinite num- 
ber of periods, i.e., as the steady state has developed. For 
large periods, the probability distribution approaches its 
equilibrium shape of a resonant-level model whereas for 
decreasing the period r side-peaks appear located at odd 
multiples of the driving frequency f2. This corresponds 
to hopping processes under the absorption or emission 
of an odd number of quanta fi. Regarding the d) and 
d operators as raising or lowering the local spin and cj_ 
as creating a spin excitation of energy k in the fermionic 
reservoir, one can deduce from Fig. [2] that the elemen- 
tary excitations caused by the periodic driving are those 
where the local spin is flipped by simultaneously creat- 
ing spin excitations of energy nQ, n odd. A Fourier series 
expansion of the #(sin(f2i))-function in the interval [0, r] 
reveals why mostly the odd frequencies contribute: 



[sin(Oi)] = - ^2 cos(nfit). 



(10) 



JME2Z+1 



Therefore, the switch on and off driving can be thought 
of as a non-monochromatic driving including all odd mul- 
tiples of the driving frequency. 

As an aside, one can deduce from Eq. Q that the 
periodic switching only leads to the enhancement of cer- 
tain fluctuations that are already present in the corre- 
sponding equilibrium Hamiltonians. The periodic time- 
dependence of the Hamiltonian leads to a selection of 
certain transitions that correspond to the absorption or 
emission of an odd number of energy quanta Q. 



IV. MAGNETIZATION 

The dynamics of all quantities is determined by the 
Hamiltonian in Eq. Q. Due to the Emery-Kivelson 
transformation, the instantaneous value of the magne- 
tization P(t) of the impurity spin 



P(t) = (S z (t)) 



(11) 



itself appears in the Hamiltonian. Therefore, the mag- 
netization provides the full access to the time evolution 
in the present time-dependent setup. As it will turn 
out in the following, this does not imply that the prob- 
lem has to be solved self-consistently. In contrast, the 
time evolution of the magnetization is independent of 
the potential scattering term appearing in the effective 
Hamiltonian as already mentioned below Eq. Q . Due to 
Eq. ([5]), the magnetization is connected to the occupation 
n d (t) = (d^(t)d(t)) of the local d-level: 



P(t) = n d (t) 



(12) 



As a consequence of the periodicity of the Hamiltonian, 
it is convenient to represent time coordinates in the fol- 
lowing way: 



t = nr + s ,n G N,s G [0,t/2]. 



(13) 



For s G [r/2, r] the S z operator is constant due to the 
switch off of the spin dynamics. Therefore, the formulas 
in the following will always be presented for s G [0, r/2]. 
Using the formulas in Eq. Q and the fact that the ma- 
trices A4^ N ^ are unitary, it can be shown that the mag- 
netization decays exponentially in time: 



P(t) = P(0) e 



-nAT„-2As 



(14) 



where the time scale is set by the Kondo scale tx — 
1/Tk — l/(wwA). For a si ngle inter action quench, the 
magnetization P qu (t) equals Bl i 6 ! 21 ! 22 ! 



P qu (t) = P(0) e 



-2At 



(15) 



Comparison with the result in Eq. ( 14 ) reveals that the 
periodic driving affects the impurity spin orientation only 
by reducing the total time during which the spin dynam- 
ics in the Kondo Hamiltonian is switched on. The initial 
local spin polarization is transferred to the conduction 
band and flows away from the central region to infin- 
ity. Note that the magnetization is independent of the 
intermediate time evolution generated by the potential 
scattering Hamiltonian. 

In the limit of fast switching, r — > 0, the magnetization 
decays exponentially in time 



P T ^o(<) =P(0) e 



-Ai 



(16) 



The associated rate A, however, is smaller compared 
to the single quench case where it is equal to 2A, see 



6 



Eq. (15 1. This is surprising, as one might expect that 



the additional energy provided to the system by the pe- 
riodic driving may open additional phase space for relax- 
ation processes. The decrease in the rate by one half oc- 
curs simply because the spin dynamics that are the only 
source of relaxation of the magnetization are switched on 
only during half of the time. 



The exact result in Eq. ( 16 ) for the magnetization, 



however, contrasts the dynamics one obtains by a naive 
implementation of the Trotter formula. Following a sug- 
gestion by Eisler and PescheP^ the dynamics of a period- 
ically quenched system in the limit r — > is identical to 
that of an effective equilibrium Hamiltonian that can be 
obtained by applying the Trotter formulsP^ to the time 
evolution operator U over one period: 



JJ — e -iH 1 T/2 e -iH 2 r/2 



■ H 1 +H 2 „ 

e - l ^~^ T + 0{t 2 ). (17) 



As a result, the effective Hamiltonian equals the time- 
averaged one, that is a resonant-level model with a hop- 
ping amplitude V/2 plus a potential scattering term. The 
potential scattering term, however, does not affect the 
spin dynamics. Therefore, it will be omitted in the fol- 
lowing. A resonant-level model with hopping amplitude 
V/2 generates a decay of the magnetization at a rate A/2 
as A a V 2 , compare Eq. (15 1, in contrast to the exact 
value A. 

From a mathematical point of view, the Trotter for- 
mula fails as both Hamiltonians Hi and H2 have to be 
self-adjoint in a mathematical sense, that is a stronger 
requirement than Hermitian. For bounded operators (in 
fact, all realistic models are equipped with a high-energy 
cutoff), self-adjointness is guaranteed, in contrast to un- 
bounded operators as considered in this work. In fact, re- 
garding a resonant-level model with a nontrivial hopping 
element V k (t) = V(t) e^l^Y where y(t) = V6[sm(Qt)} 
and k c a high-energy cutoff, one indeed observes that in 
the limit r — > the magnetization converges to the result 
obtained by applying the Trotter formula, see Eq. (17). 



Physically speaking, the inapplicability of the Trotter 
formula is a consequence of the creation of high-energy 
excitations in the fast driving limit. The typical excita- 
tions generated by the periodic driving are of the order 
of Q = 2tt/t corresponding to the absorption and emis- 
sion of quanta of the driving frequency, compare Fig. 
As £1 — > 00 for t — > 0, the typical excitations carry high 
energies. If the system is not provided with a mecha- 
nism, such as a finite bandwidth, that suppresses or cuts 
off these high-energy excitations, they will exist in the 
system even for r — > 0. Clearly, such excitations are 
not present in equilibrium systems leading to the con- 



clusion that this fast driven system cannot described by 
an equilibrium Hamiltonian reflecting the mathematical 
statement above. 

Summing up, the Trotter formula is not applicable in 
this model with a flat hybridization function. Despite 
the fact that realistic models actually exhibit a bounded 
spectrum, the results presented in this work describe cor- 
rectly the dynamics in the limit Q 3> TV, provided that 
the physical cutoff D is still much bigger than ft. 



V. SPIN-SPIN CORRELATION FUNCTION 

A dynamical quantity that carries more information 
about the local properties of the Kondo model is the 
spin-spin correlation function: 

(S z (t)S z (t'))=C(t,t')- i - X (t,t'), (18) 

where C{t,t') = \{{S z (i), S z (t')}) denotes the sym- 
metrized part and x(M') = - t')([S z (t), S z (f )]) the 
response function for t > t'. After an infinite number of 
periods r, a steady state develops in which all real-time 
correlation functions are invariant under a discrete time 
shift r in all their time arguments. For the spin-spin 
correlation function this implies 



(S z (t + r)S z (t' + t)) = (S z (t)S z (t r )). 



(19) 



Therefore, the time coordinate t' can be restricted to the 
interval [0, r]. The average (...) in the steady state is to 
be understood as 

(S z (t)S z (t')) = lim (V (NT)\S z (t)S z (t')\* (NT)), 

iv— >-oo 

(20) 

where |^o) denotes the initial state. In the steady state 
where (S z (t)) — (S z (t')) — 0, the spin-spin correlation 
function equals 



(S z (t)S z (t')) = (h d {t)n d (t')) 



1 



(21) 



Using the formulas in Eq. j9|, the spin-spin correlation 
function, that is a four-point function in terms of the 
fermionic operators of the effective Hamiltonian, can be 
related to a two-point function: 



(S z (t)S z (t')) = (S(t)d{t')) 2 . : 



(22) 



where 
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(S(t)d(t')) = 



cku 



-e- As I e 



l+UJ 2 

-iuAt' 



,-A(s+t' 



1 - 2cos(wAr/2)e- AT / 2 H 
1-2 cos(wAr)e- Ar / 2 + e~ A 



-+( 



') (< 



jAt' 



-iwA-r/2 _ -At/2 
_At'\ e e 



,-At/2 



-At' / iwAs -As' 

e e — e 



e iwAr/2 _ g-Ar/2" 



o-At/2 



(23) 



Remarkably, the real part of the (d' (t)d(t')} correlator 
can be calculated analytically 



0.25 



<{rft(t),d(f)}) 



-ijAt/2 — A(s— t') 



(24) 



matching precisely the result of Langreth and 
Nordlander 2 ^ who derived a general formula of this 
correlator for an arbitrarily time-dependent resonant- 
level model. Inserting the time-dependence of the 
hopping amplitude V(t) of the present setup into their 
result shows perfect agreement. As the anticommutator 
of two time-evolved fermionic single-particle operators is 
independent of the state if the Hamiltonian is quadratic, 
the full information about the influence of the steady 
state onto the spin-spin correlation function is contained 
in the imaginary part of the (eft (t)d(t')) correlator. The 
imaginary part is not accessible analytically, but can be 
evaluated numerically. A representative plot of C(t,t') 
and x(M') is shown in Fig. [3} A detailed discussion will 
be given below. 

Due to the appearance of a new scale in the present 
periodic driving setup, the universal description of the 
spin-spin correlation function gets modified in compari- 
son to its equilibrium form: 



(S z (t)S z (t>)) 



t t' T 
tK ' tK tK 



(25) 



as one can expect by dimensional analysis. As a conse- 
quence, the Kondo scale, though an equilibrium quan- 
tity, remains the only relevant time and energy scale as 
it was shown for the conductance thr ough a Kondo im- 
purity in the work by Kaminski et al , ! n l 12 l The period r 
only appears in combination with the Kondo time scale 
tK = defining the ratio: 



V 



T 
tK 



A" • 



(26) 



Therefore, it is only important how fast the system is 
driven compared to the internal time scale tK- 

The limit of long switching times: In the limit of large 
periods r, the behavior of the local correlation functions 
is accessible by general arguments. Roughly speaking, 
the system is able to relax during each half period. Ini- 
tially, the system is prepared in one of the ground states 
|^o) of the noninteracting Hamiltonian that are product 
states of the Fermi sea for the conduction-band electrons 
and a spin wave function for the local level. Switching 
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FIG. 3: Universal curves for the symmetrized correlation 
function C(t) = C(t, 0) and the response function \(t) = 
x(t, 0) for different values of the driving rate 77 at zero wait- 
ing time t' — 0. The parameter 77 = r/tfc compares the speed 
of the external driving r with the internal Kondo time scale 
t K = l/T K . 



on the interaction in the Kondo model creates local ex- 
citations in the vicinity of the local spin and the Kondo 
singlet forms. The excitations generated by the quench 
delocalize and flow away from the central region to infin- 
ity such that they cannot contribute to local properties 
any more. Therefore, the time evolved state looks like 
the ground state for local obscrvablcs after a sufficiently 
long time. As was emphasized in RefP, however, the 
state l^o) can never develop into the true ground state 
of the Kondo model as the overlap of both wave func- 



tions is constant in time. Nevertheless, the time evolved 
state |^o) is essentially equivalent to the true ground 
state as far as expectation values of local observables 
such as the S z operator are concerned. All statements 
about relaxation of the state itself in the following are 
to be understood in this sense. After half of the pe- 
riod, the interaction in the Kondo model is switched off, 
thereby destroying the Kondo singlet. The excitations 
that are created by breaking up the Kondo singlet delo- 
calize as argued before and the system evolves into the 
ground state of the noninteracting Hamiltonian, that is a 
product state of the Fermi sea and a spin wave function 
with zero magnetization. Therefore, the system at the 
moment of the second switch on of the interaction is pre- 
pared as initially up to a change in the local spin wave 
function. As a result, the system behaves as for a sin- 
gle interaction quench, a situation t hat has already been 
addressed by Lobaskin and Kehrein!^ Analytically, the 
spin-spin correlation function transforms into: 



{S z {t)S z {t')) 



2 



< - - i (s(t - t') 



-e~ At s(t) + e- At s{t')) , (27) 



where s(t) = ir^ 1 doj sin (a; At) /(uj 2 + 1). This result 
matches precisely the result obtained by Lobaskin and 
Kehreirp' for a single interaction quench. 

The limit of fast switching: In the opposite limit 
t — > 0, there exists no general argument capturing the 
dynamics in the resulting steady state. The sole reason 
for this is that the Trotter formula is not applicable in 
this model as explained in Sec. IV. Otherwise, the dy- 
namics in the steady state would be governed by an ef- 
fective equilibrium Hamiltonian. Nevertheless, the fast 
driving generates dynamics similar to equilibrium as the 
system is not able to follow the fast externally prescribed 
perturbation. 

Performing the limit r — t 0, the spin-spin correlation 
function reduces to: 



(s z (t)s z (t')) - 



I e -A(t-t')/2 _ i l ; 



■{{t-t')/2) 



(28) 

Thus, it only depends on the time difference signaling 
the similarity to an equilibrium problem in the sense 
that time-translational invariance is restored. The equi- 
librium spin-spin correlation function equals!^ 



(S z (t)S z (t')) eq = 



i e - A (*-*') -is(t-t') 



(29) 



Comparing the result of Eq. ( 28 ) with the equilibrium 



case, one first observes that the time argument is scaled 
by a factor of 1/2. As for the magnetization, this can 
be understood by the fact that the spin operators evolve 
nontrivially only during half of the time. Additionally, a 
prefactor of 1/2 appears in front of the function s such 



that the real and imaginary part of the (eft (t)d{t')) corre- 
lator transform qualitatively different in the limit r — > 0. 
Thus, the fluctuation-dissipation theorem, a signature of 
equilibrium systems, is violated. This then leads to the 
conclusion that it is impossible to find an equilibrium 
Hamiltonian generating the same dynamics. 

Nevertheless, other quantities such as the dynamical 
spin susceptibility, as it will be shown in Sec. VI, can be 
related to equilibrium Hamiltonians. The inapplicability 
of the Trotter formula does not exclude the possibility to 
find equilibrium behavior, it only excludes the possibility 
to find a unique equilibrium Hamiltonian describing the 
correct dynamics for all observables. Concluding, the 
steady state in the limit r — > provides an example 
for the emergence of new quantum states similar to 
equilibrium states but with new properties that are not 
accessible by equilibrium thermodynamics. 

The asymptotic long-time behavior: In equilibrium, the 
spin-spin correlation function exhibits a characteristic al- 
gebraic long-time behavior at zero temperature: 



(S z (t)S z (t')} 



t—t'- 



vq 



t-t' 



(30) 



At finite temperatures, the decay is exponential due to 
the smearing of the Fermi surface. Therefore, the long- 
time behavior of the spin-spin correlation function can 
serve as a measure whether the excitations that are cre- 
ated by the periodic driving are equivalent to those in- 
duced by a finite temperature. In this case, the steady 
state may be characterized by relating it to an equilib- 
rium system at an effective te mpera ture, a concept that 
has been widely used recentlypEEEl 

In the present periodic time-dependent setup, the 
asymptotic long-time behavior of the spin-spin correla- 
tion function can be determined analytically: 



t°o 2 
> — W 



t 



K 



t-t' 



(31) 



Surprisingly, the long-time behavior of the spin-spin cor- 
relation function is universal in the sense that it is com- 
pletely independent of the external driving. Moreover, it 
precisely matches the equilibrium behavior at zero tem- 
perature. As the algebraic decay in equilibrium is caused 
by the sharp Fermi surface at zero temperature, one can 
conclude that the periodic driving is not able to smear 
the Fermi surface or to at least locally heat up the system 
excluding the concept of effective temperature. Further- 
more, the low-energy excitations in the immediate vicin- 
ity of the Fermi level that are relevant for the long-time 
behavior are unaffected by the periodic driving. The ex- 
citation spectrum involves excitations of multiples of the 
driving frequency that emerge from processes where the 
local spin is flipped by simultaneously creating collective 
spin excitations in the fermionic reservoir with energies 
of odd multiples of the driving frequency f2, as can be 
seen in Fig. [2| Therefore, the excitation spectrum is of 
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discrete character, in contrast to the excitation spectrum 
induced by temperature. 



VI. DYNAMICAL SPIN SUSCEPTIBILITY 

In the steady state, the magnetization of the impurity 
spin vanishes due to Eq. ( 14 ) . A local spin polarization 



can be induced by applying a magnetic field to the local 
spin. In the linear response regime one obtains: 



(S,(t))h = (s z (t)) 



dt' x (t,t')h(t') 



(32) 



Expectation values without an index h are to be evalu- 
ated with respect to the unperturbed Hamiltonian. In 
the steady state, the magnetization vanishes such that 
the expectation value for the local spin polarization is 
solely determined by the response function 



X (t,t') = i9(t-t')([S z (t),S z (t')]). 



(33) 



In equilibrium, the response function x(AO only de- 
pends on the time difference thereby establishing a spec- 
tral representation of only one frequency argument whose 
imaginary part x"( e ) is called the dynamical spin suscep- 
tibility. It shows a peak located near the Kondo temper- 
ature Tk that can be associated with the Kondo singlet. 
The existence of such a peak in a non-equilibrium set- 
ting can also be interpreted as a signature for the pres- 
ence of the Kondo effect as a whole, as the Kondo sin- 
glet is just one manifestation of this coherent many-body 
phenomenon. Moreover, in equilibrium, the fluctuation- 
dissipation theorem holds that connects the dynamical 
spin susceptibility with the local spin fluctuation spec- 
trum. In systems out of equilibrium, the Fluctuation- 
Dissipation theorem is violated, as can be seen explic- 
itly for a single interaction quench scenario in the Kondo 
mo del !^ 

For periodic time-dependent Hamiltonians there also 
exists a preferable spectral decomposition. Due to the 
periodicity property in Eq. (19), it is suitable to define 
two new time arguments: 



_ t + 1' _ , 



(34) 



Expressing x(i,i') in these coordinates, x(i av ,i ro i) is in- 
variant under the transformation i av — > < av + r such that 
a Fourier series expansion in the coordinate i av can be 
performed. Therefore, one can spectrally decompose x 
in the following way that is usually referred to as the 
Wigner representation: 



Xn(e) = 

i r 



^ 1 e i£treI xKv + % i ,i a v- 



i r 



dt m e q 



2 ' 2 



i r 



dCe ,nnt "x(i„,S) 



(35) 



where e = e + nil/2. With each component n, one can 
associate the behavior of the quantity x due to processes 
where n quanta of the driving frequency il are absorbed 
(n > 0) or emitted (n < 0). The quantity x(*av> e) can be 
interpreted as the spectral decomposition of x at a given 
point i av in time. The n = component of Xn simply is 
the time average of the quantity x(^ av i s). 

Applying a small sinusoidal magnetic field h(t) = 
ho sin(fiot) to the local spin in the steady state, linear 
response theory predicts for the magnetization of the im- 
purity spin (S z (t)}h in the presence of the small pertur- 
bation h: 

(S z (t)) h = h Y, [x'ni^o + nil/2) sin[(fi + nil)t] 

n 

-Xni^o + nil/2) cos[(Q + nil)t]} . (36) 

Remarkably, this expression contains a static contribu- 
tion if ilo = mil, m € Z\{0}. Therefore, it is possible to 
align the local spin on average by applying a sinusoidal 
magnetic field that itself contains no static contribution: 



(S,(t)) h = x;(mn/2). 



(37) 



Here, (...) denotes the time average. As the spin dy- 
namics are switched off during the second half period, 
the external magnetic field is not able to influence the lo- 
cal spin magnetization in this time window regardless of 
its time dependence. Therefore, the magnetic field must 
be capable of polarizing the local spin during the first 
half period in order to induce a static component. This 
can be achieved by h(t) — ho sin(i7ot) with Qq = mil. 

Note that a definite order of taking limits is implicitly 
prescribed in this linear response calculation. First, the 
system is evolved into the steady state that is established 
after an infinite number of periods. Afterwards, an ad- 
ditional sinusoidal magnetic field is applied that acts for 
an eventually infinite amount of time. Therefore, linear 
response theory only provides the information of how 
an additional infinitesimal magnetic field influences the 
local spin properties after the steady state has already 
been established. It does not necessarily describe the 
properties of a system whose interaction is periodically 
switched on and off in presence of a magnetic field, 
as the two involved limiting processes may not commute. 



Results for x"(tm> £ ) : In Fig. El results for the dynam- 
ical spin susceptibility x"(t^,£) are shown for different 
values of the driving frequency il. The first plot dis- 
plays the case of slow driving, i.e., r>(j{, and reveals 
again the single quench dynamics as argued before. After 
a transient regime on a time scale tx, that is the time 
scale for the buildup of the Kondo effectP, the dynam- 
ical spin susceptibility approaches its equilibrium shape 
with a peak located at e ~ Tk representing the buildup 
of the Kondo singlet and therefore of the Kondo effect. 
Near the half period boundary, however, a new struc- 
ture emerges. At £ av ~ t/2 — tjc, a whole sequence of 
strong peaks appears that cannot be explained by the 
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FIG. 4: False color plots of the dynamical susceptibility 
x"(tav,e) at a given point i av in time for different values of 
the parameter r\. The upper plot shows the behavior for slow 
driving consistent with a single interaction quench picture. 
The lower plot displays the fast driving behavior that is sim- 
ilar to an equilibrium problem. The plot in the middle shows 
the intermediate regime of competing time scales. 



simple single interaction quench picture as this coher- 
ent phenomenon is solely caused by the periodic driving. 
Approaching i av = t these peaks rapidly decrease and 
disappear. Note, that the dynamical spin susceptibility 
is not discontinuous at £ av = = t/2 as it might seem. The 
collapse near t/2, however, is so fast that it cannot be 
resolved by the pictures presented here. 

In the case of fast switching, the lower plot in Fig. |4j 
the dynamical spin susceptibility stays nearly constant 



over the whole period. The system is not able to adapt 
to the fast external perturbation, as it is faster than the 
internal time scale t^, on which the system is able to 
react. Remarkably, the shape of the spectral decompo- 
sition at any time point t av resembles the shape of an 
equilibrium dynamical spin susceptibility with a reduced 
Kondo temperature Tk/2. It was shown in Eq. (28) that 
the response function Xt->-o (tj ^V) in the limit r — > 
is identical to half of an equilibrium response function 
Xeq(t/2, Tk)/2 at a rescaled time argument t/2 and iden- 
tical Kondo temperature Tr. Here, the Kondo temper- 
ature has been included explicitly for convenience. Such 
a relation for functions in time implies the following re- 
lation for the corresponding Fourier transforms: 



Xr^o(e,T K ) = Xeq{2e,T K ). 



(38) 



For the dynamical spin susceptibility, this statement can 
be rewritten in terms of a modified Kondo temperature: 

xUoi^T K ) = 2x': q (s,T K /2). (39) 

The equilibrium dynamical susce ptibil ity at zero temper- 
ature is known exactly, such tha t 16 ! 26 ! 



A 2 



1 



2tt e 2 + A 2 




1 (2e 
+ — arctan — 
A V A 



(40) 



A Kondo singlet is present even for fast periodic driv- 
ing with a reduced binding energy Tk/2. Therefore, one 
can speculate that the Kondo effect itself survives with a 
rescaled Kondo temperature Tjc / 2 Although this anal- 
ogy suggests the equivalence to an effective equilibrium 
problem, this is not valid in a strict sense, since it was 
shown, compare Eq. (28), that one cannot find an ef- 



fective equilibrium Hamiltonian that generates the same 
dynamics. Nevertheless, it is remarkable that the dynam- 
ical spin susceptibility, that is a measurable quantity in 
principle, shows an equilibrium-like behavior of a Kondo 
model with reduced Kondo temperature in this limit. 

The robustness of the Kondo effect even for ft ^> Tk 
has also been observed by Nordlander and coworkers 3 
where the time-averaged local spectral density shows a 
Kondo resonance in the limit £1 3> Tk- 

The plot in the middle of Fig. [4] shows the interme- 
diate regime where both the external time scale r and 
the internal time scale tx are of the same order. In 
this case, the system is able to adapt partially to the 
externally prescribed perturbation. For small times 
t av , one observes that the Kondo singlet tries to form, 
the dynamical spin susceptibility tends to approach its 
equilibrium profile. Near the half period boundary where 
the interaction is switched off, however, the dynamical 
spin susceptibility collapses onto a curve with a peak lo- 
cated approximately at Tk /2 as in the fast driving case. 
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FIG. 5: Universal curves of the zeroth, first and second mode 
Xn(£ — nil/ 2) for different values of the parameter n. 



Results for x^(e): In Fig. [5j universal curves for the 
modes x" (e) are shown. In the upper plot of the n = 
component, one can nicely see the crossover from the 
fast to the slow driving regime. For fast driving, small 77, 
one observes a peak located at e ~ Tr/1, whose height 
and position decreases by increasing r\. This signals a de- 



crease in the binding energy as well as the stability of the 
associated Kondo singlet. For large 77, this peak vanishes 
according to the single quench dynamics for t = 00 where 
only a Kondo singlet of energy Tk exists. Indeed, for in- 
termediate 77, a second peak emerges at Tk with a height 
of only half of its equilibrium value, as half of the time 
the interaction is switched off. Thus, the n = compo- 
nent of the dynamical susceptibility shows the crossover 
from a Kondo singlet of binding energy Tk at t = 00 to a 
Kondo singlet with binding energy Tr-/2 for r = 0. The 
oscillations appearing for large 77 originate in the periodic 
structure near the half period boundary of x(^av,£) that 
can be seen in Fig. [1] 

The n = 1 and n = 2 modes that correspond to pro- 
cesses where one or two energy quanta fl are absorbed 
are shown in the lower two plots of Fig. [5] Due to the 
relation 



(41) 



these results can be extended straightforwardly to the 
71 = — 1, —2 modes. 

For small periods r, all modes vanish due to the 

equilibrium-like dynamics where absorption and emission 
of photons is not taking place. As £1 —> 00 for r — > 0, 
the absorption or emission of an energy quantum f2 in- 
volves high energy transfers that are cut off by the finite 
width of the local hybridized d level. The equilibrium 
spin fluctuations of the conduction-band electrons in the 
Kondo model at large energies are suppressed due to en- 
ergy conservation. As mentioned in Sec. Ill, the only 
influence of the periodic driving is the enhancement of 
certain fluctuations that are present in the equilibrium 
versions of the involved Hamiltonians. If the equilibrium 
fluctuations already vanish, their enhancement due to the 
periodic driving is negligible, too. 

The larger the parameter 77, the larger becomes the re- 
sponse for processes under absorption and emission, as 
the excitations generated by the periodic driving now 
turn out to be excitations of lower energy as f2 decreases. 
In the limit r — > 00, however, the n 7^ components 
are expected to vanish as in the r — > limit, since the 
system is able to relax during each half period and lo- 
cally behaves as in equilibrium nearly all of the time. 
The behavior for the n = 1 and n — 2 components is 
qualitatively different for energies above or below nil/2. 
The absorption or emission of one photon is likely for the 
71 = 1 mode whereas the n = 2 mode decays to zero much 
faster away from its maximum. The different properties 
of the even and odd modes can be traced back to Eq. ( 10 1 , 



stating that the square wave driving can be thought of as 
a non-monochromatic driving including all frequencies of 
odd multiples of il. 



VII. CONCLUSIONS 

In this work, the steady state in a periodically driven 
Kondo model in the Toulouse limit has been investigated 
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by analyzing exact analytical results for the local spin 
dynamics at zero temperature. Remarkably, the anal- 
ysis revealed a universal asymptotic long-time behavior 
of the spin-spin correlation function, compare Eq. (311, 



independent of the driving matching precisely the zero 
temperature equilibrium result. This universality origi- 
nates in the discrete character of the excitation spectrum 
that leaves the low-energy excitations in the immediate 
vicinity of the Fermi level unchanged. 

For large switching times, r/tx — ► oo where is the 
equilibrium Kondo time scale, the local observables be- 
have as for a single interaction quench, since the system 
is able to relax during each half period. In the opposite 
limit, r/tx — > 0, the system is not able to follow the 
fast externally prescribed perturbation. Although the 
dynamics is "equilibrium-like" (in the sense that time- 
translational invariance is restored, e.g. Eq. (28)), it is 



impossible to find an effective equilibrium Hamiltonian 
generating the same dynamics. This can be seen in the 
discussion below Eq. (29) were it was shown that the 



fluctuation-dissipation theorem is violated for the local 
spin observable. Consequently, the steady state created 
by the periodic time-dependent setup provides an exam- 
ple for the possibility to reach new quantum states in sys- 
tems out of equilibrium that are not accessible by equi- 
librium thermodynamics. The reason for the existence 
of this unusual steady state is the inapplicability of the 
Trotter formula in the present periodic time-dependent 
setup. This observation might also be useful for other 
periodically driven systems. 

Remarkably, in the Toulouse limit the Kondo effect is 
robust against the periodic driving as long as the driving 
frequency is unable to induce charge fluctuations in the 
quantum dot. Even in the case of fast driving Q 3> Tjf 
the Kondo singlet survives. The associated binding en- 
ergy Tx/2 is reduced to half of its equilibrium value. For 
intermediate r, a crossover is observed from a Kondo sin- 
glet of binding energy T^/2 in the fast switching limit to 
a Kondo singlet of binding energy Tk for long switching 
times. 

It would be interesting to investigate related questions 
for the Kondo model away from the Toulouse limit, where 
one also expects that the flow to the strong-coupling fixed 
point itself is affected by the periodic switching. Likewise 
it would be important to understand whether other pe- 



riodically driven quantum impurity models show similar 
dynamics. Work along these lines is in progress. The 
exact analytical results derived in this work will be im- 
portant benchmarks for these investigations. 
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Appendix A 

In this appendix, the procedure to obtain the solution 
to the matrix multiplication problem is presented. For 
that purpose, the Green's functions for the resonant-level 
model Q and the potential scattering Hamiltonian C/W 
have to be derived, for example, by using the equations 
of motion approach. The matrices Ai(n) — QQ^ n > obey 
the following formulas: 



M d d{n) = 
Mdk{n) = 

M kd (n) = 



V 
k-iA 

y e -ikr/2 



p -ikr/2 _ e -Ar/2 



e -ifer/2 _ e -Ar/2 



k - iA 

Mkk' (n) = 8kk'e~ lkT + £ kk ' + K>kk> (n) 
where the matrices C and IC(n) are determined by: 

e -ifcr/2 



(Al) 



y2 e -ikr/2 



(k - k')(k - iA) 



-At/2 



-ik'r/2 



/C fc fc'(n) 



(k - iA)(k' - iA) (k' - k)(k' - iA) 
g(S z (nT))e- tkT / 2 e~ lkT ' 2 - e r ik ' T / 2 



l + ig(S z (nT))L/2 



k-k' 



(A2) 



For simplicity, the method to obtain will be dis- 

played by calculating the matrix element M dd • Using 
its definition in Eq. (fSl) one can write: 



M$ = ]T M dh (l)...M lN _ ld {N)=M ( d N d 1] M dd + Yl M dll (l)...M lN _ 3d (N-2)M dkN _ 2 M kN _ 2d 
+ E Mdi 1 W---M lN _ 3 k N _ 2 (N-2)M kN „ 2 k N „ 1 (N-l)Mk N _ 1 d. (A3) 

I 

To proceed further, the following relations are required: that originate in elementary properties of the matrix M, 

namely that A4 as a function of k exhibits no poles and 

Y,M k y (n)M k , d e- iXk ' r = e-^ +1 ^M kd , 

k' 

J2M dk M kd = (A4) 

k 
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that the exponentials e~ lkT / 2 always allow to deform in- 
tegration contours into the lower half plane. Using these 
formulas, one can directly show that the two sums ap- 
pearing in Eq. ( A3 1 vanish. The resulting recursion for- 



mula is solved by: 



M 



dd 



-NAt/2 



(A5) 



Analogously, the other matrix elements of can be 

derived, although the evaluation of M^J is quite in- 
volved. Remarkably, all matrix elements with an index 



d are independent of the potential scatterer. Only M^O 
includes a term that can be traced back to the potential 
scattering: 



™ 1 „—inkT~—i(N—ri)k'T 



n=0 



■ig(S,(0)W 



(A6) 



Here, the result for the magnetization, see Eq. (14 1, has 
been used. 



1 B. P. Anderson and M. A. Kasevich, Science 282, 1686 
(1998); M. Greiner, O. Mandel, T. W. Hansen and I. Bloch, 
Nature 419, 51 (2002). 

2 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, 
D. Abusch-Magder, U. Meirav and M. A. Kastner, Nature 
391, 156 (1998); S. M. Cronenwett, T. H. Oosterkamp and 
L. P. Kouwenhoven, Science 281, 540 (1998); J. Schmid, 
J. Weis, K. Eberl and K. von Klitzing, Physica B 258, 
182 (1998): W. G. van der Wiel, S. De Franceschi, T. Fuji- 
sawa, J. M. Elzerman, S. Tarucha and L. P. Kouwenhoven, 
Science 289, 2105 (2000). 

3 P. Nordlander, N. S. Wingreen, Y. Meir and D. C. Lan- 
greth, Phys. Rev. B 61, 2146 (2000). 

4 D. Lobaskin and S. Kehrein, Phys. Rev. B 71, 193303 
(2005). 

5 D. Lobaskin and S. Kehrein, J. Stat. Phys. 123, 301 (2006). 

6 M. Moeckel and S. Kehrein, Phys. Rev. Lett. 100, 175702 
(2008); M. Eckstein, M. Kollar and P. Werner, Phys. Rev. 
Lett. 103, 056403 (2009); P. Calabrese and J. Cardy, Phys. 
Rev. Lett. 96, 136801 (2006); C. Kollath and A. M. Lfiuchli 
and E. Altman, Phys. Rev. Lett. 98, 180601 (2007); 
A. Faribault, P. Calabrese and J.-S. Caux , J. Stat. Mech., 
P03018 (2009); M. Eckstein and M. Kollar, Phys. Rev. 
Lett. 100, 120404 (2008); M. Rigol, Phys. Rev. Lett. 103, 
100403 (2009); S. R. Manmana, S. Wessel, R. M. Noack 
and A. Muramatsu, Phys. Rev. Lett. 98, 210405 (2007); 
M. A. Cazalilla, Phys. Rev. Lett. 97, 156403 (2006). 

7 P. Nordlander, M. Pustilnik, Y. Meir, N. S. Wingreen, 
D. C. Langreth, Phys. Rev. Lett. 83, 808 (1999). 

8 A. Mitra and A. J. Millis, Phys. Rev. B 76, 085342 (2007). 

9 A. V. Joura, J. K. Freericks and Th. Pruschke, Phys. Rev. 
Lett. 101, 196401 (2008); N. Tsuji, T. Oka and H. Aoki, 
Phys. Rev. Lett. 103, 047403 (2009); N. Tsuji, T. Oka and 
H. Aoki, Phys. Rev. B 78, 235124 (2008). 

R. Lopez, R. Aguado, G. Platero, C. Tejedor, Phys. Rev. 
Lett. 81, 4688 (1998); T. K. Ng, Phys. Rev. Lett. 76, 487 



(1996); M. H. Hettler and H. Schoeller, Phys. Rev. Lett. 
74, 4907 (1995); B. H. Wu and J. C. Cao, Phys. Rev. B 
77, 233307 (2008). 

A. Kaminski, Y. V. Nazarov and L. I. Glazman, Phys. Rev. 
Lett. 83, 384 (1999). 

A. Kaminski, Y. V. Nazarov and L. I. Glazman, Phys. Rev. 
B 62, 8154 (2000). 

Y. Goldin and Y. Avishai, Phys. Rev. Lett. 81, 5394 
(1998); Y. Goldin and Y. Avishai, Phys. Rev. B 61, 16750 
(2000). 

A. Schiller and S. Hershfield, Phys. Rev. Lett. 77, 1821 
(1996). 

G. Toulouse, C. R. Acad. Sci. Paris 268, 1200 (1969). 

A. J. Leggett, S. Chakravarty, A. T. Dorsey, 
M. P. A. Fisher, A. Garg and W. Zwerger, Rev. 
Mod. Phys. 59, 1 (1987). 

B. Doyon and N. Andrei, Phys. Rev. B 73, 245326 (2006). 
J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 
(1966). 

J. von Delft and H. Schoeller, Ann. Phy. (Leipzig) 7, 225 
(1998). 

G. Zarand and J. von Delft, Phys. Rev. B 61, 6918 (2000). 
F. Lesage and H. Saleur, Phys. Rev. Lett. 80, 4370 (1998). 
F. Guinea, V. Hakim and A. Muramatsu, Phys. Rev. B 
32, 4410 (1985). 

V. Eisler and I. Peschel, Ann. Phys. (Berlin) 17, 410 
(2008). 

H. F. Trotter, Proc. Amer. Math. Soc 10, 545 (1959); 
M. Suzuki, Commun. Math. Phys. 51, 183 (1976). 

D. C. Langreth and P. Nordlander, Phys. Rev. B 43, 2541 
(1991). 

F. Guinea, Phys. Rev. B 32, 4486 (1985). 
Note that the time averaged effective Hamiltonian would 
generate a Kondo model with a Kondo temperature Tk /4 
in contrast to the exact value Tk/2. 



